Surface deformation caused by the Abrikosov vortex lattice 
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In superconductors penetrated by Abrikosov vortices the magnetic pressure and the inhomoge- 
neous condensate density induce a deformation of the ionic lattice. We calculate how this deforma- 
tion corrugates the surface of a semi-infinite sample. The effect of the surface dipole is included. 
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I. INTRODUCTION 

Deformations of the ionic lattice caused by the 
Abrikosov vortices have been studied from several as- 
pects. For example, if a vortex moves, the motion of such 
deformation demands a motion of ions which contributes 
to the inertial mass of the vortex.i^ In the static case, 
forces evoked by vortices create a tension which modifies 
the total volume of the sample and which was observed 
as magnetostriction.— Moreover, in anisotropic materials 
the elastic energy caused by vortices depends on the rel- 
ative orientation of the Abrikosov vortex lattice and the 
crystal lattice. 5 ' 6 

In the above mentioned studies a vortex was treated as 
infinitely long. This idealized geometry essentially sim- 
plifies the problem. Since the system has translational 
symmetry along the vortex line, the lattice deformation 
is purely longitudinal with the displacement vectors per- 
pendicular to the vortex* 1 ^ 

As far as we know, nobody has studied deformations 
near the surface, where the magnetic flux of the vortex 
leaves the superconductor. In the present paper we fo- 
cus on this problem. For simplicity we assume that the 
sample is semi-infinite and the applied magnetic field is 
perpendicular to its surface, see figure [T] 

Far from the surface the lattice deformation ap- 
proaches its bulk value with displacement vectors per- 
pendicular to vortices, i.e., parallel to the surface. Let 
us sketch in advance which complications one can expect 
near the surface. 

As seen in figure [TJ the stretching magnetic field to- 
gether with circulating currents result in a Lorentz force 
with a component pointing along the vortex. Electric 
fields balancing this force act also on the ion lattice and 
corrugate the surface. 

Beside the Lorentz force, one can imagine two addi- 
tional mechanisms leading to the corrugation. First, the 
specific volume of the superconducting and normal states 
differ. The normal metal in the vortex core is pulled by 
its neighborhood to adopt the specific volume of the su- 
perconductor. As mentioned above, deep in the bulk 
these stresses cause displacements perpendicular to the 
vortex. Close to the surface, however, the same stresses 
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FIG. 1: Lorentz force acting on the circulating superconduct- 
ing current. In the bulk the Lorentz force is parallel to the 
surface. Near the surface it is not parallel due to the magnetic 
stray field. 



lead to different deformations as they are partially re- 
laxed by displacements perpendicular to the surface, i.e., 
along the vortex. 

Second, the surface dipole which confines the electrons 
in the crystal changes under transition from the normal 
to the superconducting state. Accordingly, at the vortex 
core the surface dipole differs from other regions of the 
surface. The force which holds electrons inward natu- 
rally pulls the oppositely charged ions outwards and also 
corrugates the surface. The general relation between the 
surface dipole and the surface tension has been estab- 
lished already two decades agoi^ The contribution of 
the superconducting condensate has been discussed only 
recently^ 

The paper is organized as follows. In sections II and 
III we introduce the basic set of equations for elastic de- 
formations of an isotropic material near the surface. The 
surface deformation due to the Abrikosov vortex lattices 
is discussed in section IV and the numerical solutions 
for Nb and YBa2Cu3C-7 are presented in section V. Sec- 
tion VI with the conclusions ends the paper. 
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II. ELASTIC DEFORMATIONS IN ISOTROPIC 
MEDIUM 

In this section we recall the basic definitions needed to 
describe the elastic deformations of an isotropic medium. 
More details the reader can find in the textbook of Lan- 
dau and Lifshitz. 10 



A. Tensors of strain and stress 

Deformations are described by atomic displacements 
u. We assume that all atoms are identical and only a 
single atom occupies the elementary cell. Since the inter- 
atomic distance is short compared to the characteristic 
scales of deformations, we treat u(x, y, z) as a function 
of continuous coordinates. 

The space derivatives of the displacement define the 
strain tensor. For small deformations assumed here its 
components read 

1 / du x du y \ 

The complementary anti-symmetric combination of 
derivatives corresponds to the rotation of the rigid body, 
therefore it does not contribute to deformations. 

The strain creates a stress described by a tensor a. In 
general a KT — Yluv ^-kt/muU^ui where A is 4 th order tensor 
of elastic coefficients. The Greek indices stand for x, y, z. 
For an isotropic medium this relation simplifies to 

a KT = K (V • u) 8 KT + (u KT - i (V • u) S KT ^j , (2) 

where the divergence of the displacement vector 

(V • u) = u xx + u yy + u zz (3) 

is a shorthand notation for the trace of the strain tensor. 
It represents the local change of the specific volume. The 
tensor u KT — |(V • u) S KT has zero trace and describes 
purely shear deformations. The coefficients K and fi are 
called the bulk and shear modulus, respectively. 



B. Stability conditions 

The stress depends on local gradients of the displace- 
ments, i.e., on changes of the bond lengths between neigh- 
boring atoms. Accordingly it represents only contact 
forces while long-range forces have to be covered sepa- 
rately. Here we shall consider a non-contact force due 
to interaction of the crystal lattice with superconducting 
electrons. 

The gradient of the stress balances a long-range force 
F acting on a unitary volume ^ K V K <r KT + F T — 0, where 



V x = gj; and so on. For the isotropic material this sta- 
bility condition reads^ 

(k + V(V • u) - M [V x [V x u]] = F. (4) 

In this paper we consider forces that can be expressed 
as a gradient of the potential U 

F = -VU. (5) 

In our numerical studies below we assume that this po- 
tential is exclusively due to the electrostatic potential (p 
acting on the charge density of the ionic lattice p, i.e., 
U = pip. 

Alternatively, one can assume an effective potential 
U = 2Ka 1 ip\ 2 jn corresponding to the force used by 
Simanck and Duani^ and Coffey^. Here oto is the relative 
volume difference between the superconducting and the 
normal states, K is the bulk modulus, ip is the Ginzburg- 
Landau (GL) wave function and n is the density of pair- 
able electrons. 

These two choices of the potential U yield very simi- 
lar results. They are not completely identical, however. 
For example, the electrostatic potential known as the 
Bernoulli potential covers the Lorentz force, while the 
effective potential U — 2Ka \t/j\ 2 jn does not^ 

C. Surface conditions 

At the surface the balance of forces demands 
~Y^ K n^KT + Pt = 0, where n is a unit vector normal 
to the surface and P is the external force on a unit area 
of the surface. The normal component p = (n ■ P) is the 
pressure. We will not assume tangential surface forces in 
this paper. 

To simplify boundary conditions we specify the geom- 
etry of our sample. It is a semi-infinite superconductor in 
the half-space z > 0. For the normal vector n = (0, 0, —1) 
and normal surface force P = (0, 0, —p) we find three sur- 
face conditions 




Using the linear relation between stress and strain, the 
surface conditions ([6]) are converted into conditions for 
the displacement u. For the isotropic medium from ([2|) 
and J6J follows 

u xz = 0, 

Uyz = 0, (7) 

(k + ^fij u zz + (k - ^f^J (u xx + Uyy) = -p. 

The pressure p includes the ambient pressure of coolant 
which is homogeneous all over the sample surface. We ig- 
nore it, because it does not contribute to the corrugation. 
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We do not want to list all possible forces acting on the 
surface. We merely mention that the surface of the sam- 
ple feels the surface dipole 8<p, which acts as an effective 
pressure p = —pSip, see Ref. d. The dipole amplitude 
Sep = cp(0) — (f(oo) is the difference between the electro- 
static potential at the surface and deep in the bulk of the 
superconductor. The surface dipole is conveniently eval- 
uated from the GL free energy with the help of the Budd- 
Vannimenus theore m 12 ' 13 adopted to superconductor^^. 

Now the problem is fully specified. The displacement u 
is driven by the potential U and the effective pressure p. 
These quantities can be treated within any approxima- 
tion the reader prefers. In section V we will relate them 
to the electrostatic potential tp and the surface dipole Sip. 
At the moment we merely assume that the potential and 
the pressure are known functions of the GL wave func- 
tion. Since the deformation has only a negligible effect 
on the GL wave function, the GL wave function obtained 
from the ordinary GL theory can be used. We take the 
GL wave function as known and focus on the deforma- 
tion. 



III. INDUCED AND FREE DEFORMATION 

As mentioned in the introduction, the deformation 
near the surface has three sources - magnetic stray field, 
expansion or contraction along the vortex, and the sur- 
face dipole. Since the elastic equation ((4]) and its bound- 
ary conditions ([7]) are linear in u, it is possible to separate 
individual contributions. 

As the first step, it is advantageous to separate the 
displacement u into two parts - induced deformation u 1 
driven by the long range force F which induces displace- 
ments all over the sample, and the free deformation u f . 
So we write 

u = u' + u f , (8) 
where the displacement u 1 is longitudinal 

[V x u 1 ] = 0. (9) 
and obeys the equation 

(k+^ (V-u i ) = -U(r,z). (10) 

Since the differential equation p0[) is of the first order, 
this induced displacement is fully specified by the poten- 
tial U and the requirement of convergence deep in the 
sample. 

After substitution of ([8]) into the elastic equation |4}, 
the stress from the induced displacement compensates 
the long range forces and one is left with the equation 
for the free displacement 

(k + ^fij V(V • u f ) - n [V x [V x u f ]] = 0. (11) 



This is a second-order differential equation and we need 
to find the solution which in addition to the convergence 
deep inside guarantees also the fulfillment of the surface 
boundary condition. 

A. Laplace equation for the free deformation 

Now we show that it is possible to simplify equation 
(fTTjl to the Laplace equation. 
Taking the divergence of equation (fTTj) we find that 

V • (V 2 u f ) = 0, (12) 

because the divergence of a rotation is zero. Taking the 
rotation of equation (|11[) we find 

[V x (V 2 u f )] =0, (13) 

because the rotation of a gradient is also zero. The vector 
identity [V x [V x u f ]] = V(V • u f ) — V 2 u f was used in 
the rearrangement. 

According to (fT2")) and (JT3J) all the derivatives of V 2 u f 
are zero, therefore it is a constant. A finite value of u f 
for z — ► oo is possible only if this constant is zero. We 
thus have 

V 2 u f = 0. (14) 

Being free of material parameters, this Laplace equation 
is more convenient than the equation pip . 

B. Surface matching 



The surface boundary condition Q applies to the total 
displacement. Substituting @ into ([7]) we obtain 




These conditions determine gradients of u at z = 0. 
The induced deformation u 1 is already known, therefore 
we have moved its strain elements to the right hand sides 
of the boundary conditions. 

By conditions (fT5|) the free displacement is fully spec- 
ified. The second boundary condition is the request of 
convergency, u l — » for z — > oo. Now the set of equations 
is complete and we are ready to solve for the corrugation 
of the surface around the vortex ends. 



IV. VORTEX LATTICE 

If the vortices form an Abrikosov lattice, we can benefit 
from its periodic structure. In this case all components 
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of the displacement can be solved analytically in the two- 
dimensional (2D) Fourier representation. 



A. Induced deformation 

In the x—y plane, the potential is a sum of planar waves 
£/(r,z) = ]Te* kr [/(k,z), (16) 

k 

where r = (x, y) is a 2D coordinate. The 2D wave vectors 
k = (k x , k y ) attain discrete values given by the density of 
vortices and the actual structure of the Abrikosov lattice. 

To solve the equation for the displacement requires 
some intermediate steps. Since the induced displacement 
u 1 is longitudinal, it can be expressed as the gradient of 
a scalar function 

u' = V X (r,4 (17) 
In the 2D Fourier representation the equation (fTO)) yields 
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k»-JL) x =U[k.z). (IN) 



where k 2 = k 2 +k 2 . It is solved by 



X(k, z) 



v 

1 1 

2k K + §/z 



U(k,z') 



(19) 



We are interested namely in values at the surface 

oo 

1 1 



X(k,0) = 



2k K + |/i 



dz'e~ kz U(k,z'). 



(20) 



The z-gradient at the surface follows from (fl"9|) as 

= fe%(k,0). (21) 



oz 



z=0 



The surface value of this displacement contributes to the 
surface corrugation. 

The surface strain given by the second gradient follows 
from the equation (fig)) as 



d 2 



fc 2 X (k,0) 



1 



2 = 



K 



-U(k,0). 



This surface strain enters the boundary condition, where 
it serves as the source of the free deformation. 



The boundary conditions (I15|) in the 2D Fourier rep- 
resentation read 



-{ik x u{ - kill) = -ik x kx{k,0), 
2 V 



-{ikyu{ - ku y ) = -ik y kx(k,0), 



(24) 



-k ( K + ^ u[ + \ K - ^m) (ik x ul + ik y u y ) 
= -p-2 A ifc 2 x(k,0) + C/(k,0). 



We have used equation (1231) to evaluate the derivatives 
on the left hand sides of (fl5|) and the ansatz (fl7| together 
with relations (|20H22p in the right hand sides of (Ti~5)) . 

From the first and second conditions of |24]) one elimi- 
nates u x and u y obtaining ik x u l x + ik y u y = —ku f z — 2k 2 \- 
Using this relation in the third condition of (f2"4"| one ar- 
rives at the z-component of the free displacement 



p-U(k,0)-2k 2 (K-| M ) X (k,0) 
2k(K+\ii) 



(25) 



With u\ from (f25|) and the first and the second equa- 
tions of (f2"4"]l one readily evaluates the parallel displace- 
ments u x and u y , respectively. Here we want to concen- 
trate on the component u l z perpendicular to the surface. 



C. Surface corrugation 

The corrugation of the surface is given by the total 
displacement in the z-direction u z = u z + u\. Adding 
formulas (|2"5|) and (f2"Tj) we obtain 



u»(k) 



« 2 r 



p - U(k, 0) + 4fc 2 /x x(k,0) 
2k (K + i/x) 

«.(k). 



E^ kr 

k 



(26) 



Formula (|2^|) is the final result of the general part of 
our discussion. It provides the 2D Fourier decomposi- 
tion of atomic displacements at the surface layer. In our 
notation a positive displacement u z corresponds to the 
atomic motion inwards the superconductor. 



V. NUMERICAL RESULTS FOR THE 
ELECTROSTATIC INTERACTION 



B. Free deformation 

The free deviation u f has the same periodicity as the 
Abrikosov vortex lattice. From the Laplace equation |TJ 
thus follows 



u'(k,z) = u t (k,0)e- fc 



kz 



(23) 



To proceed we assume that the interaction between the 
super-electrons and the ionic lattice is exclusively carried 
by the mean electrostatic field. This approximation ne- 
glects the effect of the lattice density on the electronic 
band structure and on the phonon band structure and for 
many materials it might lead to quantitatively incorrect 
predictions. In this paper we adopt this approximation 
for its simplicity. 
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A. Surface dipole 

As mentioned in the introduction, in the electrostatic 
approximation, the effective potential is proportional to 
the electrostatic Bernoulli potential tp, i.e., 

U(t,z) = (hp(t,z). (27) 
The Bernoulli potential is a function of the GL function 15 



etp 



da 
dn 



10/3 
20^ 



M 4 



(28) 



where a = 7 (T 2 - T c 2 )/2n and = jT 2 /2n 2 are the GL 
parameters expressed in terms of Sommerfeld's gamma. 
These values follow from the Gorter-Casimir two-fluid 
model in the limit T — > T c , see Ref. 16. One can 
also use the limiting BCS values of a and (3 which are 
larger by factors 1.43 and 2.86 than the two-fluid limit, 
respectively.— 

The pressure at the surface is caused by the surface 
dipole^ 



p(r) = p<p+(v) -p<p_(r), 



(29) 



where <p±(r) = ip(r, ±e) with e > being an infinitesimal 
distance. In reality, e has to exceed the Thomas-Fermi 
screening length and the BCS coherence length. The 
potential ip+ is the extrapolation of the internal potential 
towards the surface, and tp_ represents the potential out 
of the crystal. 

The surface value of the effective potential is defined 
as a limit from inside, therefore 



U(r,0) 



limZ7(r,e) = pp+(r). 

e~ >0 



(30) 



The displacement of the surface atoms (|26l) for the elec- 
trostatic approximation reads 



2k (K 



2kp 



K 



- X (k,0). (31) 



As one can see, the potential ip + cancels and we are 
left with the potential tp_. According to the Budd- 
Vannimenus theorem for superconductors, the potential 
above the surface is proportional to the density of the 
free energ y 14 i 17 



e<f- 



2n 



(32) 



which one conveniently evaluates within the Ginzburg- 
Landau theory. For details of the Bernoulli potential and 
the surface dipole see textbook [l6| . 



B. Bulk forces 

The internal value of the Bernoulli potential con- 
tributes only via the subsidiary function \. In general 
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TABLE I: Material parameters of Niobium and YBa2Cua07: 
condensation energy per particle, GL parameter, particle den- 
sity and logarithmic derivatives of the critical temperature 
and the linear coefficient of the specific heat 7, Young's mod- 
ulus and Poisson ratio. 



one needs the complete 3D solution of the vortex lattice 
near the surface. Numerical studies have shown that one 
can neglect the z dependence of the amplitude of the GL 
function everywhere including the close vicinity of the 
surface. 18 The amplitude variation is about 1%. Since 
the Bernoulli potential depends exclusively on the am- 
plitude of the GL function, the potential has the same 
value as deep in the bulk, tp(r, z) ~ (^^(r), for any z. 

Neglecting the z dependence of the potential, the inte- 
gral in (j2"0|) becomes trivial giving 



X(k,0) 



1 



P 



2k 2 K 



¥>oo(k). 



(33) 



The potential ip^ is conveniently found from (|28p using 
the GL function deep in the sample. For a fixed magnetic 
field, the GL function deep in the sample is the same as 
in the infinite sample, therefore it can be evaluated as an 
effectively two-dimensional problem. 



C. Surface corrugation 

The displacement of the surface atoms results from 
3D and (ESS) as 



2k (K + \p) 



p_(k) 



2ft 



K 



4 <Poo(k) 
3M 



(34) 



One can see from (|34|) that the relative contribution of 
the bulk and the surface potentials depends on the shear 
modulus p. For soft materials with p <g; K the bulk po- 
tential gives a negligible contribution so that the surface 
corrugation is driven only by the electrostatic potential 
above the surface. In this case the corrugation is inde- 
pendent of the density derivatives of the GL parameters. 

Many authors dealing with deformable superconduc- 
tors prefer to express elastic coefficients via the Poisson 
ratio <7 = [K — / (2A^ + and the Young modu- 



lus E = 3K(1 - 2a). 
of surface atoms (l34l) 



In this notation the displacement 
reads 



(1-2<t)(1+<t) 
kE(l-a) 



{(l-(r)p(k)-^oo(k)}. (35) 



Deriving (|35j) we have used that the Bernoulli potential is 
nearly independent of z so that <p + — ip^. The pressure 
at the surface (|29[) thus equals p = p{poo — <f— ) ■ 
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FIG. 2: Surface corrugation u z (r) of Nb (above) compared 
with the value neglecting surface dipoles (middle) and only 
surface dipoles (below). The temperature and the mean mag- 
netic field are T = 0.95 T c = 9 K and B = 0.21 B c2 {T) = 
6.4 mT. Length unit is the vortex distance a = 128 ran. 

The numerical study of the surface corrugation (|35p is 
performed with the parameters listed in the table Q] Fig- 
ure presents the surface corrugation of Niobium with 
the GL parameters k increased by impurities to 1.5. The 
scale is in pm=10 -2 A, therefore the maximum deforma- 
tion ~ 5 10~ 4 A is far too small to be detected by recent 
scanning microscopes. In our convention the negative it z 
corresponds to atoms displaced out of the crystal. 

The total deviation of atoms at the surface is seen in 



Uzl 



u z [nm; 



u z [nm; 




y -o.5 



FIG. 3: Surface corrugation u z (r) of YBa2Cus07 (above) 
compared with the value neglecting surface dipoles (mid- 
dle) and only surface dipoles (below). The temperature and 
the mean magnetic field are T = 0.67 T c = 60 K and 
B = 0.01 B c2 {T) = 0.6 T. Length unit is the vortex distance 
a — 58 nm. 



the upper part of figured The middle part shows the de- 
viation evaluated omitting the surface dipole, i.e., setting 
p = in formula (|35[) . As one can see, such approxima- 
tion for Niobium leads to the opposite sign of the atomic 
displacement. The lower part shows the surface deforma- 
tion evaluated from the surface dipole only. This approx- 
imation overestimates the amplitude more than twice. 
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The magnitude of the surface corrugation is directly 
proportional to the condensation energy per particle, 

which for Niobium has the value = 4.6 /icV. In 
conventional superconductors this condensation energy 
per particle is small because of the small critical tem- 
perature, T c = 9.5 K and the large density of particles, 
n = 2.2 10 28 /m 3 . 

In high-T c superconductors the critical temperature is 
larger by a factor of ten, while the density of particles 
is lower by a factor of ten. This leads to an appreciably 
larger condensation energy per particle. For example, in 
YBa 2 Cu 3 7 one has T c = 90 K and n = 5 X 10 27 m" 3 
resulting into = 750 fieV. Figure [3] presents a corru- 
gation of the YBa2Cu307 surface with the surface par- 
allel to the a-b planes, i.e., the magnetic field along the 
c axis. As one can see, the maximum displacement of 
the surface atom is more than 0.5 A, which can be ob- 
served by the scanning microscope. Again, neglecting the 
surface dipole one arrives at the opposite displacements, 
while using only the surface dipole one overestimates the 
magnitude. 

Perhaps we should note that application of the 
isotropic model to the layered structure of YBa 2 Cu307 is 
not justified. One clearly needs more elastic coefficients 
to describe the deformation of this highly anisotropic 
material^ The presented data can thus serve only as 
an order of magnitude estimate. 



VI. CONCLUSIONS 

Magnetic field entering the superconductor in the form 
of vortices induces a corrugation of the surface. In con- 
ventional superconductors the displacement of surface 
atoms is of the order of 10~ 4 A, which is too small to be 
observed by recent experimental tools. In high-Tc super- 
conductors one can expect amplitudes ~ 10 _1 A which 
is in the reach of scanning force microscopes. 

Our results for Niobium and YBa2Cu307 show that 
among the forces that drive the surface corrugation the 
dominant one is due to the surface dipole. The contri- 
bution of the bulk potential to the surface corrugation is 
opposite to the contribution of the surface, therefore it 
reduces the magnitude of the atomic displacement. 

Since the deformation of the crystal near the surface 
differs from the deformation in the bulk, one can expect 
that the surface terms play an important role in the 
Simanek contribution to the vortex mass in thin layers. 
We leave this problem for a future work. 
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